WHAT IS CLAIMED IS: 



1 . A method of magnetic gradient wavefomi design for rapid magnetic resonance 
imaging using constrained optimization comprising the steps of: 

5 a) selecting a niunber N of discrete-time waveforms with a sampling period , x, and 

time T=Nt, 

b) expressing constraints as equations, 

c) identifying if a solution exists for T, and if not, increasing T until a solution exists, 

d) decreasing T to find a shortest solution, and 
10 e) solving the equations to identify the solution. 

2. The method as defined by claim 1 wherein step e) utilizes Unear programming to 
identify the solution. 

15 3. The method as defined by claim 2 wherein the linear programming finds the vector x 
that minimizes a cost function f ^x subject to the constraint Ax = < b, where matrix A and 
vector b are formed by combining all of the linear constraint equations for the amplifier and 
pulse sequence constraints. 

20 4. The method as defined by claim 1 wherein steps c), d), and e) utilize LI -norm 
formulation. 

5. The method as defined by claim 4 wherein the gradient is one dimensional. 
25 6. The method as defined by claim 4 wherein the gradient is two dimensional. 

7. The method as defined by claim 4 wherein the gradient is three dimensional. 

8. The method as defined by claim 4 wherein in step a) the gradient amplifier has current 
30 and voltage limits (Imax and Vmax) that result in the following limits on G(p): 

\G (01 < ri W 

and 
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L~-G{t)-^RG{t) 
at 



< vV 



where L, R, and r| are the gradient coil inductance, resistance, and efficiency. 



9. The method as defined by claim 1 wherein steps c), d), and e) include adding slack 
5 variables to the optimization. 

10. The method as defined by claim 9 wherein step e) solves for gradient Gx(n) and a set 
of variables Hx(n) that converge to |Gx(n)| and a set of variables Sx(n) that converge to 
|aGx(n) + /3Gx(n+l)|. 

10 

1 1 . The method as defined by claim 10 wherein the slack variables are forced to converge 
by adding linear gradient constraints: 

^//,M-fG,w<o 

and slew-rates constraints: 

j3 " SArt] + aG^n] + >6GJ« + 1] < 0 

- S^n] - aG,[n] - pG^n + 1] < 0 

whereby the constraints force the Hx(n) and Sx(n) variables to approach an 
appropriate absolute value when combined with minimization of the following pulse 
function: 

/I=l 

20 

12. The method as defined by claim 1 wherein step e) utilizes second-order cone 
programming (SOCP) to find a solution x that minimizes a linear pulse fiinction f ^x subject 
to a second order cone constraint 

II^jc + ZjII^ <Cx + d 

25 which is a superset of the linear constraints in linear programming where matrix A 

and vector b are formed by combining all linear constraint equations for amplifier and pulse 
sequence constraints. 
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13. The method as defined by claim 1 wherein steps c), d) and e) use a binary-search to 
minimize time in identifying the minimum value of n which provides a solution. 



14. The method as defined by claim 13 wherein the gradient is one dimensional. 

15. The method as defined by claim 13 wherein the gradient is two dimensional. 

16. The method as defined by claim 13 wherein the gradient is three dimensional. 

1 7. The method as defined by claim 13 wherein in step a) the gradient amplifier has 
current and voltage limits (W and V^ax) that result in the following limits on G(p): 

|G(0|<r|Ux 



and 



L — G{t) + RG{t) 
dt 



— 'I max 

where L, R, and r| are the gradient coil inductance, resistance, and efficiency. 

18. The method as defined by claim 1 wherein the gradient is one dimensional. 

19. The method as defined by claim 1 wherein the gradient is two dimensional. 

20. The method as defined by claim 1 wherein the gradient is three dimensional. 

21. The method as defined by claim 1 wherein in step a) the gradient amplifier has current 
and voltage limits (Imax and V^ax) that result in the following limits on G(p): 

|G (/)|<r|Ux 

and 

j^G(t) + RGit)<r,V^ 
where L, R, and r| are the gradient coil inductance, resistance, and efficiency. 



22. The method as defined by claim 1 and fUrther including: 
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f) establishing physical constraints of a gradient amplifier and a gradient coil. 

23. The method as defined by claim 22 and further including: 

g) establishing beginning and end gradient boundary constraints and fixed area for 
one or more gradient moments. 

24. The method as defined by claim 1 and further including: 

f) establishing beginning and end gradient boundary constraints and fixed area for 
one or more gradient moments. 

25. The method as defined by claim 1 wherein step b) includes gradient heating 
constraints. 

26. The method as defined by claim 25 wherein step b) includes gradient 
magnetostimulation constraints. 

27. The method as defined by claim 1 wherein step b) includes gradient 
magnetostimulation constraints. 
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